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INTRODUCTION 

Much  of  the  recent  work  on  the  statistical  mechanics  of  the  electrical  double  layer  has 
focused  on  the  application  of  theories  based  upon  integral  equations  to  models  where  the 
solvent  is  assumed  to  be  structureless.  These  theories  are  an  improvement  of  the  classical 
Poisson-Boltzmann  theory  of  Gouy  [1]  and  Chapman  [2],  mainly  because  they  give  a  better 
treatment  of  the  short-ranged  part  of  the  ionic  interactions.  However,  the  influence  of  the 
solvent  structure  on  the  electrode-electrolyte  interface  is  well  recognized  [3].  It  is  thus 
desirable  to  extend  the  integral  equation  approach  to  models  where  the  solvent  is  another 
molecular  species. 

'  In  this  communication -we  consider  a  simplified  model  of  an  electrode-electrolyte 
interface  which  takes  into  account  the  discreteness  of  the  solvent  (Fig.  1).  The  electrode  is 
represented  as  a  hard  infinite  plane  with  a  uniform  surface  charg?  density,  while  the 
electrolytic  solution  is  modeled  as  a  mixture  of  charged  hard  spheres  (the  ions)  and  dipolar 
(nonpolarizable)  hard  spheres  (the  solvent).  In  addition,  we  assume  the  electrode  has  a 
dielectric  constant  unit;  thus  we  ignore  the  image  effects.  For  this  model  we  developed  a 
treatment  in  which  all  the  interactions  (ion-ion,  ion-solvent,  solvent-solvent,  electrode 
solution)  are  considered  on  the  same  (classical)  statistical  mechanical  basis.  In  particular  the 
solvent  is  taken  as  a  fully  orientable  system  and  no  special  treatment  of  a  hypothetical  inner 
layer  is  made. 

i  ■ 

Blum  and  Henderson  [4]  and  also  Camie  and  Chan  [S]  have  studied  this  model  within 
the  mean  spherical  approximation  (MSA),  the  simplest  known  integral  equation.  We  have 
considered  elsewhere  [6]  the  application  of  another  integral  equation,  related  to  the  MSA, 
which  may  be  referred  as  a  generalized  MSA  (GMSA).  In  this  approximation  the 
Ornstein-Zemike  (O-Z)  equation  for  the  particle-particle  and  electrode-particle  correlations 
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is  solved  for  closures  which  are  generalizations  of  those  of  the  MSA.  The  correlation 
functions  thus  obtained  give  information  of  the  micro-strucrure  of  the  system  and  also  permit 
the  evaluation  of  several  thermodynamic  and  electrostatic  properties  of  the  interface.  Here 
we  present  some  further  GMSA  calculations  of  some  of  these  properties  which  complement 
those  of  Reference  [6]. 

THEORY 

a.  Electrode-Particle  Correlation  Functions 

We  divide  the  total  (h)  and  direct  (c)  electrode-ion  and  electrode-dipole  correlation 
functions  into  a  nonelectrostatic  and  an  electrostatic  part 

h±(x)  -  h^il^x) 

hd(x,0)  -  h^x)  +  /3  hd(x)  cos  0  (la) 

c±(x)  -  Cq(x)  ±Cj(x) 

cd(x,0)  -  c0(x)  +  vT  cd(x)  cos  9  ,  (lb) 

x  and  9  being,  respectively,  the  reduced  distance  (x/D)  from  the  electrode  and  the  angle 
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and  we  assume 

c0(x)  -  A0  e“ V 
Cj(x)  »  E*q*x  +  ^-*1*  , 

cd(x)  -  +  Ad  e-V  ,  x>i  ,  (3) 

where  Aa,  za  (a-S4,d)  are  parameters  to  be  determined.  When  As-Ai«Ad-0  we  recover 
the  MSA  closures  [4].  In  Eq.  (3)  we  use  the  dimensionless  quantities  q*^»/Jq^/D; 

and  E*»(/JD^)1/2e  (E«4iro),  /3  being  the  Boltzmann  thermal  factor  (kBT)“*. 

If  we  consider  row  matrices  Mw[M|,Md],  then  the  O-Z  equation  for  the  wall-particle 
correlation  functions  becomes  [6,7] 

b(x)  -  f(x)  +  f*  dt  b(x  -  t)QB(t)  ,  x>i  ,  (4) 

where 

1  (x)  *  dk  e_ik*  £(k)  [5b(  -  k)  ]  (5) 

with 

C(k)  -  /J2  dx  e^  £(x)  (6) 

and 

QbW -I-JJ'dre^  SB(r).  (7) 

Here  QB(r)  is  the  Baxter  matrix  of  the  electrolyte  and  is  of  the  form 

fQu(r)  Qid(r)l 

SbW  - 

LQdi(r)  Qdd(f)J 

Similar  relations  (involving  scalars)  are  valid  for  the  nonelectrostatic  terms  [8]. 


(8) 


We  have  obtained  explicit  expressions  for  Qg(r)  in  the  MSA  [9,10]  and  GMSA  [11]  as 
simple  sums  involving  powers  of  r  and  exponential  functions  of  r.  If  we  use  the  Eq.  (3)  with 


Aa  -  0  (a— s4,d)  as  the  closure  of  the  O-Z  equation  [Eq.  (4)]  we  have  a  MSA/MSA  or  a 
MSA/ GMSA  calculation  according  to  whether  we  use  the  bulk  MSA  or  GMSA  expressions 
for  Qb. 

If  we  consider  instead  Aa,  za  (a«s,i,d)  as  adjustable  parameters  we  have  a 
GMSA/MSA  calculation  when  Qg  (MSA)  is  used  or  a  GMSA/GMSA  calculation  if  Qg 
(GMSA)  is  utilized.  In  both  these  last  cases  we  can  obtain  the  parameters  demanding  that 
the  correlation  function  satisfy  the  linearized  hypemetted  chain  (LHNC)  [12]  conditions 


h«(x“  i)  “expl 

.’>.(*- i)]-1  (<»-+.-) 

(9a) 

M>  - 1) - 

(9b) 

>>'.(*  -  -  1.(x  ■  j)«p[>I.(*  -  1)]  (a  -  +  ,-) 

(9c) 

-  \  )  -  -  -j)  +  [",(*  -  y)  +  >K(X  -  t)  • 

(9d) 

where 

ly(x)  -  hr(x)-cT(x)-/3UT(x)  (y  -  +  ,-,d) 

(10) 

with 

0U+ 

(x)  -  -£U_(x)  -  -E*q*x 

(11a) 

P*M* 

W"  -  -7T 

(lib) 

and 

&(*-  i)  ■ 

-  t(*  - 1) 

(12) 

h'(x  -  T )  -  t  (*  "  2  )  +  "  2  )Qb  (r  “  0> 

(13) 
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k.  Smm  Straetaral  and  Electrostatic  FncdoM  of  tho  Interface 

Once  we  obtain  the  total  wall-particle  correlation  functions  by  numerical  integration  of 
Eq.  (4)  (adequate  precision  is  obtained  using  the  trapezoidal  rule),  we  can  evaluate 
1  —  The  density  profiles  of  the  ions  and  the  density-orientation  profile  of  the  dipoles 


Pt( x)  -  p±[h±(x)  +  1] 
Pd(x,0)  -  pd[hd(x,0)  +  1]  , 


(14) 


2  -  The  charge  density  profile 


3  - 


q*(x)  -  p* q*hj(x)  , 
where  p*  «  pjD3  , 

The  polarization  density  profile 


(15) 


P*(x)  -  (0D)1/2P(x)  -  hd(x)  , 

v  3 


(16) 


4  - 


with  pd  »  pdD3,  and 

The  potential  profile  4  *(x)«(/fD)*/2$(x),  the  solution  of  the  Poisson  equation. 


(17) 


which  is  given  by 


**(x)  -  — 4irp*q*J^"  dt  (t  -  xlh^t)  -  ~£r  pd p*  f*  dt  hd(t) 


A**-E*x  +  4wp,*q*J^*  dt  (t  -  x)h,(t)  +  pdp*f^  dt  hd(t)  ,  (18) 


where  A$*-<p*(x*0)  is  the  total  potential  drop  across  the  interface. 


!  1 


j 


er> 


We  obtain  [6]  analytical  expression*  for  A4*  and  also  for  P*  the  total  polarization 


P*  -  03D)1/2P  -  /’  dt  hd(t) 

v3  J0 


Hm  hd(s)  , 

S-+0 


&<t>*  -  -4rp,* q* J^*dt  thj(t)  -  -£=•  pJm*  dt  hd( 


-4»  M»)]. 


,(s)  -  — 


[a^id  +  a^O-Tij)] 


Sd^  dQfl-Tdd.!) 


2  ,  ,  n\  '  { [al<w] ,  >1-  tsJo*.)]  =4 1 . 

(■i,d  ^«dg(l  *  ddjj)  1  ad  ai  J 


(s)  E 


*■*0  as  , 


dj  <2 


,  [d  “Tji)(l  —  Tdd)— TjrfTdi] 


flod-W 


+  «Sd  VM1 


^r— r  j  [fil(j*.)]  (1  "  Tdd)  +  [5b 

-  tdd;l^  t  «« 


and  j  *  v-1  and  at  and  aj  are  the  non-null  elements  of  the  matrix 


Pi  a*l 

^  a  1  ■  -  lim  Qt 

Lo  (J  r*- 


and  Tjj.  Tdd.  etc.,  are  the  elements  of  the  matrix 


I  -  P  dr  [QB(r)  +  AJ  ■ 


... ; 
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The  bulk  strength  parameters  and  d2  are  defined  by 

dj  -  4irp*q*2  ;  d|  -  -y-  pJm*2  • 


Finally 


(l  -  TW  -  [(1  -  T^)2  +  (T^-df] 


From  Eq.  (20)  we  can  evaluate  the  differential  capacity  from 


_ 

d  d  4ir  4v'  dE*  ' 


RESULTS 

In  Figs.  2  and  3  we  show  the  charge  and  polarization  density  profiles  calculated  by 
using  the  GMSA/GMSA,  GMSA/MSA  and  MSA/MSA  approximations.  The  oscillations 
observed  are  consequence  of  the  discrete  character  of  the  solvent.  This  follows  from  a 
simple  observation  of  the  corresponding  curves  resulting  from  models  where  the  solvent  is 
treated  as  a  continuum  [13,14]. 


In  Fig.  4  we  compare  the  potential  profile  for  the  three  integral  equations  considered. 
Here  again  the  oscillations  are  notable.  The  change  in  slope  of  the  potential  at  the 
Helmholtz  plane  (x- 1/2),  which  implies  a  change  of  the  direction  of  the  electric  field,  is  due 
to  the  large  polarization  at  this  plane.  This  is  more  clearly  seen  if  we  write 


d»*(x) 

dx 


E*  — 4irP*(x  -  I)  . 


If  we  attempted  to  define  a  local  dielectric  constant  in  the  impenetrable  region  or  Stern 
layer  [0sxil/2]  by  requiring  that  the  electric  displacement  be  continuous  at  x-1/2  we 
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would  arrive  at  a  negative  local  dielectric  constant  in  this  region.  This  illustrates  the 
difficulties  and  ambiguities  associated  with  attempts  to  define  a  local  dielectric  constant.  The 
point  is  that  the  dielectric  constant  is  a  macroscopic  quantity,  computed  using  an  infinitely 
large  volume,  and  any  attempts  to  define  a  local  dielectric  constant  will  lead  to 
inconsistencies. 

Figures  5  and  6  show  the  differential  capacity  as  a  function  of  the  total  potential  drop 
across  the  complete  double  layer.  The  MSA/MSA  curve  is  independent  of  A<£  and  so  is 
unphysicai  except  near  A<£— 0.  The  other  two  curves  are  more  reasonable  and,  in  particular, 
show  the  differential  capacitance  increasing  with  the  magnitude  of  the  potential.  Ultimately, 
we  expect  GMSA  differential  capacitance  curves  to  flatten. 

CONCLUSIONS 

The  MSA  gives  considerable  insight  into  the  structure  and  properties  of  the  double 
layer.  However,  it  has  the  drawback  that  the  MSA  results  are  linear  in  the  charge  density  on 
the  electrode.  The  GMSA  results  presented  here  overcome  this  difficulty.  These  results 
represent  our  first  attempt  to  apply  the  GMSA  and  have  some  difficulties.  For  example,  the 
density  profiles  do  not  satisfy  the  contact  value  conditions  of  Blum  and  Henderson  [4].  We 
expect  to  develop  other  and  better  schemes  for  choosing  the  GMSA  parameters  in  the  near 
future. 

In  contrast  to  the  featureless  curves  of  the  Gouy-Chapman  theory,  the  MSA  and 
GMSA  ionic  profiles,  potential  profile  and  polarization  density  show  oscillations  which  imply 
a  layering  of  sheets  of  alternating  charge  near  the  electrode.  These  oscillations  are  even 
more  pronounced  in  the  GMSA  than  in  the  MSA.  In  contrast  the  differential  capacitance 
curves  are  rather  featureless.  In  view  of  the  simplicity  of  our  model  which  ignores  the 
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interaction  of  the  metal  with  the  double  layer  and  which  is  based  oq  a  oversimplified  picture 
of  the  solvent,  the  differential  capacitance  curves  are  reasonable,  especially  considering  that 
we  have  not  attempted  to  adjust  any  parameters  to  fit  any  experimental  data. 
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